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Method and arrangement for reducing a number of 
measured values of a technical system 

The method relates to a method and an arrangement for 
reducing a number of measured values of a technical 
system. 



10 A technical system, for example a process engineering 
plant, supplies a large quantity of measured values per 
unit time, using different measured-value pickups 
(sensors) . In the course of a number of days or weeks, 
a quantity of data accumulates which requires a 

15 correspondingly high computational power to process it. 
If the measured values are to be used to adapt or reset 
the technical system, operations are often necessary 
whose complexity admits only some of the measured 
values. It is then a great disadvantage to select a 

20 specific number of measured values arbitrarily from the 
quantity of all measured values and to process them 
further, since measured values with low significance, 
for example measured values with a high degree of 
interference, have a considerable influence on the 

25 overall result and falsify the latter considerably. 



The object of the invention is to specify a method and 
an arrangement for reducing measured values, it being 
ensured that the remaining ' measured values have a high 
30 significance with regard to their description of the 
technical system. 



35 



This object is achieved in accordance with the features 
of the independent patent claims. Developments of the 
invention emerge from the dependent claims. 
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in order to achieve the object, a method of reducing 
number of measured values of a technical system 
specified, in 
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which the measured values are divided into classes in 
accordance with predefined criteria. The measured 
values in a class are assessed and those measured 
values whose assessment lies below a predefined first 
5 threshold value are screened out. 

Screening out the measured values results in a 
reduction with regard to the number of measured values.. 
A considerably reduced number of measured values is 
10 thus available for further processing. The further 
processing can be carried out with a lower 
computational outlay as compared with the non-reduced 
number of measured values. 

™ 15 One development consists in that the classes themselves 

are also assessed. In particular, a class whose 

-'M assessment lies below a predefined second threshold 

IP 

value can be screened out. This results in an 
^ additional reduction in the number of measured values . 

ffl 20 

I Another development consists in that one criterion for 

G the division into classes consists in that, for each 

class, measured values relating to a predef init ion of 
■== setting parameters of the technical system are 

O 25 determined. Typically, the technical system is set 

^ using a predefined number of setting parameters and 

after the setting has been carried out, there follows a 
(mostly time-delayed) reaction of the system to the 
setting parameters (transient response, transient 
30 process of the system) . After the setting operation, a 
specific quantity of measured values which can be 
associated with the transient process are therefore 
picked up, measured values which are associated with 
the predefined set of setting parameters continuing to 
35 accumulate after the transient process has concluded 
(transition to steady-state operation) . By adjusting 
the setting parameters, a new class is determined. All 
the measured values which respectively accumulate after 
an 
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adjustment to the setting parameters belong in their 
own class. 
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Another development consists in that measured values 
within a class which can be associated with the 
respective transient process are screened out. In 
addition, erroneous measured values can be screened 
out. In many cases, the setting of large technical 
systems depends on long-term, steady-state operation. 
Measured values which relate to the transient process 
(of short duration in comparison with the steady-state 
operation after the transient process has concluded) 
are expediently screened out, since measured values for 
the . steady-state operation are falsified by them. In 
particular within the context of modeling the technical 
system, the measured data about the steady-state 
behavior of the technical system are of interest. 

One refinement consists in reducing the number of 
measured values in a class in that at least one 
representative value for the measured values in the 
class is determined. Such a representative value may 
be: 

a) an average (e.g. a sliding average) of the 
measured values of the class, 

b) a maximum value of the measured values of the 
class , 

c) a minimum value of the measured values of the 
class , 

d) a median. 

In the case of variant d) , one advantage resides in the 
fact that it is always possible to determine a value 
which is actually present, while the average a) does 
not itself occur as a value. 

Depending on the application, a suitable selection for 
determining the representative value of a class may be 
made . 



GR 98 P 5866 

- 3a - 

An entire class with measured values can be screened 
out if said class contains less than a predefined 
number of measured values. 

Another refinement consists in that those measured 
values are screened out which differ from a 
predefinable value by more than a predefined threshold 
value . 
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The predefinable value can be an average of all the 
measured values of the class or a measured value to be 
expected in response to the respective setting 
parameters of the technical system. 

A' development consists in that the measured values 
which are reduced in number are used for the simulation 
and/or the draft design of this or another technical 
system. 

The draft design of the technical system can in this 
case comprise both adaptation and redrafting of this or 
another technical system. 




15 In addition, within the context of the simulation, the 
behavior of the technical system can be checked by 
using the reduced number of measured values, with the 
objective of altering the system for developing a new 
system with altered predef initions . 

20 

In addition, in order to achieve the object, an 
arrangement for reducing a number of measured values of 
a technical system is specified which has a processor 
unit, which processor unit is set up in such a way that 
25 the measured values can be divided into classes in 
accordance with predefined criteria. Measured values in 
a class can be assessed, measured values whose 
assessment lies below a predefined first threshold 
value are screened out. 



30 



This arrangement is particularly suitable for carrying 
out the method according to the invention or one of its 
developments explained above. 



35 



Exemplary embodiments of the invention will be 
illustrated and explained below using the drawing, in 
which 
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Fig. 1 shows a block diagram which contains steps in 
a method for reducing measured values; 

Fig. 2 shows a schematic sketch of a recovery 
5 boiler; 

Figs. 3-5 show input variables, actuating variables and 
outward variables of the recovery boiler. 

10 Fig. 1 shows a block diagram which contains steps in a 
method of reducing measured values. In a first step 
101, the measured values are divided into classes. In 
particular, the division into classes is carried out 
with the effect that all the measured values which 

15 belong to one set of setting parameters are combined in 
one class. To this extent, each alteration to the 
setting parameters of the technical system found a new 
class. In particular, altering the setting parameters 
is associated with a transient process of the technical 

20 system, this process, as opposed to a steady-state 
behavior, containing extreme fluctuations in the 
measured values. In a step 102, individual measured 
values within one class are screened out. These may be, 
for example, erroneous measured values, that is to say 

25 measured values which exhibit a high deviation with 
respect to the other measured values or an average of 
the measured values or measured values from the 
transient process. There are a number of possible ways 
of screening out individual measured values within a 

30 class: 

1. Measured values which are too poor (based on a 
predefined comparative value) ; 

2 . Measured values which belong to a transient 
process ; 

35 3. Measured values which are erroneous; 

4 . Determining a representative measured value as 
a representative for a plurality of measured 
values, in that the representative measured 
value is determined as an average of the 
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s or as a maximum 
n^easured values in a c as measured 
value or a minxmum value 

values . 
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Measured values of this type are preferably not taken 
into account; they are removed from the respective 
class. This results in a considerable reduction in the 
number of measured values. In a step 103, individual 
5 classes are screened out. One criterion for screening 
out an entire class consists in that the class contains 
less than a predefined number of measured values. In a 
step 104, the measured values reduced in number are 
used for further processing. Further processing is, in 
10 particular, a simulation and/or a draft design of the 
technical system. 



Fig. 2 shows a schematic sketch of a recovery boiler. 
In the following text, by using .the example of a 
15 ^'recovery boiler", an exemplary embodiment of the 
method described above will be illustrated. 

^1 In the paper and pulp industry, various chemicals and 

^ also heat and electrical power are needed for the 

^ 20 digestion of pulp. With the aid of the recovery boiler, 

the chemicals used and additional thermal energy may be 
Q recovered from a thickened waste process liquor (black 

^= liquor) . A measure of the recovery of the chemicals is 

of critical importance for the economy of the overall 
D 2 5 plant. 

The black liquor is burned in a char bed 201. In the 
process, an alkaline melt is produced and flows away 
via a line 202. In further process steps, the chemicals 

30 used are recovered from the constituents of the 
alkaline melt. Heat of combustion which is released is 
used to generate steam. The combustion of the waste 
liquor and therefore the recovery of the chemicals 
begins with the atomization of the black liquor via 

35 atomizer nozzles 204 into a combustion chamber 203. As 
they fall through the hot flue gas, particles of the 
atomized black liquor are dried. The dried liquor 
particles fall onto the char bed 201, first combustion 
and 



40 
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chemical reduction taking place. Volatile constituents 
and reaction products pass into an oxidation zone, in 
which oxidizing reactions proceed and in which the 
combustion is completed. 

5 

Important objectives for the control of the recovery 
boiler are the steam production in order to obtain 
power, compliance with emission values from 
environmental points of view and the efficiency of the 
10 chemical reduction. 

The combustion operation, and therefore the objectives, 
are controlled in particular by the supply of air at 
three levels (primary air (PA) , secondary air (SA) , 
15 tertiary air (TA) ) . The overall process is subject to 
numerous influences, which have to be taken into 
account during the modeling: 

a) the measurement of the variables is subject to 
20 fluctuations which are often extreme; 

b) influencing variables which are not measured 
and cannot be measured exist; 

c) at each alteration to the settable parameters 
transient processes occur; 

25 d) the technical plant becomes soiled and is 

cleaned at predefined intervals, which has the 

effect of a drift over time in each case in the 
system behavior. 

30 The measured variables of the overall process are 
subdivided into input variables (cf. Fig. 3) and output 
variables (cf. Fig. 5). Measured variables are stored 
every minute. Four of the input variables are 
simultaneously also actuating variables (also: settable 

35 parameters; cf. Fig. 4). The actuating variables are to 
be viewed substantially as free parameters of the 
overall process which can be set independently of one 
another. Some of the other 
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input variables are more or less dependent on the 
actuating variables. According to one predef inition, 
the variables "BL Front Pressure" and "BL Back 
Pressure'' are always to be regulated equally in the 
5 recovery boiler. The four actuating variables (cf. Fig. 
4) are preferably to be stored as actuating variables 
(with the desired, preset value) and as input variables 
(with the measured, actual value) . 

10 In the recovery boiler, one problem consists in the 
fact that, depending on the settable parameters, 
specific objectives, which are defined via measured 
variables, have to be met. Here, a three-stage 
procedure is selected in order to solve the problem: 

The objectives to be considered are modeled by 
means of stochastic methods, these models being 
updated by means of new measurements (data-driven, 
empirical modeling) . In this case, it is expedient 
to use not just a single model but global models 
for the identification of interesting areas in a 
parameter space determined by the objectives, and 
to use local models for the exact calculation of 
optimum operating points. The models used are 
assessed by means of quality measures . 

If the models considered are not sufficiently 
accurate because of the state of the data (quality 
measure) , new operating points are deliberately 
evaluated in order to improve the model 
(experimental design) . In addition, by using 
global stochastic optimization methods with regard 
to the objectives, attractive regions are 
identified on the basis of the current global 
model . 

For the local optimization, local models are 
constructed, and the data sets which are available 
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are, if appropriate, deliberately expanded 
(experimental design) , 
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The objectives are physical/technical or economic 
criteria which, as a rule, have to meet boundary 
conditions and/or safety conditions. It is often the 
case that a number of these criteria have to be 
5 considered at the same time. A stochastic model can be 
used in particular for the purpose of simulating the 
objectives to be optimized and their dependence on the 
parameters to be set in the computer. This is necessary 
when measurements are very costly or very time- 
10 consuming. In the case of safety requirements, possible 
hazardous situations can be avoided. 

In the case of the recovery boiler, on-line 
optimization, which is based on a plurality of items of 

15 data, is necessary, since the physical/chemical 
processes cannot be modeled quantitatively with 
sufficient accuracy and because the behavior of the 
plant is subject to fluctuations in the course of 
operation. The knowledge about this behavior must 

20 continually be expanded by means of the deliberate 
selection of new operating points. Therefore, within 
the context of on-line optimization, the above- 
described three-stage procedure of stochastic modeling 
at mathematical optimization is to be recommended. 

25 

DESCRIPTION OF THE INPUT VARIABLES 

The a input variables (a e N, N: set of natural 
numbers) generally depend on n actuating variables n e 
30 N and on random effects. They can be described as 
follows : 

Let (Q, 5, P) be a probability space and be a Borel 
a-algebra over (R: set of real numbers) for each v g 
35 N. The input variables are represented by a projection 
(p which can be measured via x S - B^: 



cp: R"" X Q -> R^ 



(1) 
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The definition set of the projection cp is a cartesian 
product of two sets. If one considers the respective 
projections onto the individual sets, then the 
following projections are obtained: 

(Pj^: n R^, CD ^ (p(x, (o) for all x € (2), 

cp^:R^ R^ X -> (p(x, co) for all co e Q (3) . 

{(Px; X e R^^} is a st ochastic process having an index set 

and a projection cp*" is a path in this stochastic 
process for each event co e Q. 

O In the case of the recovery boiler, n=4 and a=14 

^ (following the elimination of the variable ^'BL Back 

Ifl Pressure") . 

J 15 Because of the required ability to measure the 

£1 projection cpx for each x e R^, the projection cpx is a 

B random variable. Under suitable additional 

preconditions, expected values and higher moments can 
be considered. This access makes the step possible from 
20 stochastic models to deterministic optimization 
problems. In the case of a deterministic optimization 
problem, the target function can be set directly by 
means of a variable, while the stochastic variable 
influences the target function but does not permit any 
25 deliberate setting. 



10 



DESCRIPTION OF THE OUTPUT VARIABLES 



The process model M of the recovery boiler will be 
30 described as a function depending on the input 
variables and further random effects. In this case, let 
(Q, Sf P) be the above probability space. The process 
model M is then a projection which can be measured by 
x S - B^: 



35 
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M : X n (4) / 

where b refers to the number of output variables. 

Since the recovery boiler is subject to a cyclic drift 
over time (from cleaning phase to cleaning phase) , a 
description using a time parameter is also conceivable. 
The output variables may be represented by projections 
v|/ that can be measured by x S - B^: 

\}/ : R^ X n -> R^ (5) , 

(x, o) M(q)(x, co), co) (6). 

If the respective projections onto the individual sets 
of the definition set are considered, then the 
following projections are obtained 

vj/x: n R^, CO ^|/(x, Q)) for all x e R^ (7), 
M/^:R^ -> R^, X M/(x, cn) for all co G Q (8) . 

|v|/x; X £ R^l is a stochastic process having an index 

set R", and the projection Vj;"* is a path in this 
stochastic process for each (O e Q. 

In the recovery boiler, b=15. 



The fact that, when defining vj/, no distinction is drawn 
between the events co used, does not mean that there is 
25 any restriction, since Q can be represented as a 
cartesian product of an Qi and an Q2 • The above 
representation therefore also comprises the model: 



10 
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(x, 0)1, 02) m((p(x, ©1), 02) (10). 
DESCRIPTION OF THE DATA SETS AVAILABLE 

Using the descriptions in the two preceding sections, 
it is possible to combine the input variables and the 
output variables together to form measured variables O. 
O is a projection that can be measured by x S - EP, 
where m = a + b, and 

CD : X O (11) , 



(x, co) 



cp(x, co)^ 



(12) 



If the respective projections onto the individual sets 
of the definition set are considered again, then the 
15 following projections are obtained: 

<I>x: n -> R"^. o> 0(x, co) for all x e R^ (13), 
O^rR^ R^, X -> ct)(x, 0)) for all co g Q (14) . 

|<t>x; X e R^l is a stochastic process with an index set BP" 

and the projection O*" is a path in this stochastic 
20 process for each co e Q. 

For each chosen tuple x of actuating variables, a large 
number of implementations of Ox in the recovery boiler 
are determined and stored, that is to say for each Xj e 
25 R", numerous implementations 
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screened out, since these two values must be equal 
according to the predef inition of the plant 
control system. The loss of data is very small. 



5 2 



10 



The data sets are divided into classes in which 
the four setting parameters (PA, SA, TA, BL Front 
Pressure, see above) are successively constant 
over time, that is to say the jth class consists 
of the data sets Dj,. 

Classes in which there are fewer than 30 data sets 
are screened out, in order that transient 
processes do not have any great influence. 



15 4. For each class, an arithmetic average Oj and an 
empirical standard deviation sj are determined for 
all the measured variables: 



1 ^ 



(17) , 



m 



20 



25 



I ^ k=l J 



(18) 



Classes in which the averages for the variables 
PA, SA, TA or BL Front Pressure are too far 
removed from the corresponding setting parameters 
are screened out. In these classes, therefore, the 
setting values could not be reached. 
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CHARACTERISTIC STATISTICAL VARIABLES FOR THE GIVEN 
CLASSES AND THEIR GRAPHICAL REPRESENTATION 



In addition to the arithmetic averages and the 
empirical standard deviations which have been 
determined for the individual classes, a common 
standard deviation s is further determined in 
accordance with 



j=i 



^2 



1 

^2 



(19) 



Here, u stands for the number of classes (205 here) and 
V for the sum of Vj, that is to say v is the number of 
all the measured values used (38,915 here). 



